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Influence of desert dust intrusions on ground-based and satellite- 
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[i] The desert dust aerosols strongly affect propagation of solar radiation through the 
atmosphere, reducing surface irradiance available for photochemistry and photosynthesis. 

This paper evaluates effects of desert dust on surface UV erythemal irradiance (UVER), as 
measured by a ground-based broadband UV radiometer and retrieved from the satellite 
Ozone Monitoring Instrument (OMI) at Granada (southern Spain) from January 2006 to 
December 20 1 0. The dust effects are characterized by the transmittance ratio of the 
measured UVER to the corresponding modeled clear sky value. The transmittance has an 
exponential dependency on aerosol optical depth (AOD), with minimum values of M3.6 
(attenuation of ~40%). The OMI UVER algorithm does not account for UV aerosol 
absorption, which results in overestimation of the ground-based UVER especially during 
dust episodes with a mean relative difference up to 40%. The application of aerosol 
absorption post-correction method reduces OMI bias up to ~13%. The results highlight 
great effect of desert dust on the surface UV irradiance in regions like southern Spain, 
where dust intrusions from Sahara region are very frequent. 
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1. Introduction 

[ 2 ] Mineral dust aerosol plays an important role in Earth’s 
climate system, absorbing solar and thermal radiation and 
modulating Earth’s radiative budget. The main sources of 
mineral dust are the desert areas, with the Sahara being the 
most important source in the northern hemisphere [Prospero 
et al. , 2002; Ginoiix et al. , 2004; Papayannis et al., 2005; Liu 
et al. , 2008]. Air masses loaded with Saharan dust particles 
frequently reach Spanish and Portuguese regions. Several 
studies analyzed Saharan dust contribution to ambient levels 
of suspended particulate matter, studying the synoptic 
meteorological conditions responsible for the transport of the 
dust air masses [Rodriguez et al. , 2001; Escudero et al, 
2005, 2006; Querol et al ., 2009]. Other studies focused on 
the retrievals of micro-physical and optical properties of 
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Saharan dust using passive remote sensing measurements with 
sun-sky photometers [Alados-Arboledas et al. , 2003, 2008; 
Lyamani et al., 2004, 2005, 2010; Elias etal., 2006; Toledano 
et al., 2007; Cachorro et al. , 2008; Prats et al., 2008; Wagner 
et al., 2009; Valenzuela et al., 2012a, 2012b]. Lidar systems 
have also been used to characterize the vertical profile and 
structure of desert dust plumes [Perez et al. , 2006; Guerrero- 
Rascado et al., 2008, 2009; Cordoba-Jabonero et al, 2011; 
Preifiler et al., 2011]. 

[ 3 ] However, there are relatively few studies analyzing 
effects of dust intrusions on shortwave solar radiation reaching 
the Earth’s surface. [Diaz et al., 2001; Lyamani et al, 2006; 
Santos et al. , 2008; Cachorro et al. , 2008; Anton et al , 2012a]. 
To our knowledge, only Diaz et al. [2007] and Anton et al. 
[2012b] have analyzed the atmospheric aerosol effects on 
spectral UV irradiance during two Saharan dust events in 
South Spain. In general, there are only a few works about this 
subject in literature [e.g., di Surra et al, 2002; Meloni et al., 
2003; Kalashnikova et al., 2007; Garcia et al, 2009] due to 
the scarcity of routinely operational ground-based stations 
with high-quality instrumentation to measure simultaneously 
UV irradiance and aerosol data during desert dust intrusions. 

[4] The analysis of the diverse atmospheric factors 
affecting the UV irradiance is motivated by the harmful 
effects of this radiation on human health, ecosystems, and 
materials [Diffey, 1991, 2004]. This paper focuses on the 
analysis of the influence of desert dust aerosol on the UV 
erythemally weighted surface irradiance (UVER) measured 
at Granada, a non-industrialized medium-sized city in 
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southeastern Spain. The study analyzes desert dust intru- 
sions detected during a period from January 2006 to 
December 2010 to evaluate the differences between the 
UVER measurements and the satellite retrievals from Ozone 
Monitoring Instrument (OMI) on NASA EOS Aura satellite 
[Tanskanen et al., 2007]. Previously, largest differences 
between the satellite and ground-based UV irradiance data 
were reported in urban polluted areas, with elevated levels of 
UV-absorbing aerosols [McKenzie et al ., 2001; Kazantzidis 
et al ., 2006; Tanskanen et al ., 2007; Ialongo et al ., 2008; 
Buchard et al. , 2008]. Thus, it is expected that desert dust 
particles (with significant absorption in the UV spectral 
range) will also produce large differences between the sat- 
ellite-derived and surface measured UV irradiance [e.g., 
Anton et al. , 2012b]. 

[ 5 ] The article has been organized as follows: Section 2 
describes the instruments and data used in this study. 
Section 3 explains methodology. Section 4 discusses results 
and section 5 summarizes main conclusions. 

2. Instruments and Data 
2.1. Ground-Based Measurements 

[6] The experimental data used in this study have been 
collected at the radiometric station located on the rooftop of 
the Andalusian Center for Environmental Studies (CEAMA, 
37.17°N, 3.6 1°W, 680 m a.s.l.) in Granada, southeastern 
Spain. This station is operated by the Atmospheric Physics 
Group (GFAT) of the Granada University. 

[7] A broadband UV radiometer, model UVB-1, manu- 
factured by Yankee Environmental Systems, Inc. (Massa- 
chusetts, U.S.), measures spectrally integrated UV irradiance 
weighted with the erythemal action spectrum adopted by the 
Commission Internationale de l’Eclairage (CIE) [McKinlay 
and Diffey, 1987] (denoted as UVER). Measurements were 
sampled every ten seconds and recorded as one minute mean 
voltages on Campbell CR10X data acquisition systems. 
Output voltages are converted into UVER values applying 
the calibration factors derived from two calibration cam- 
paigns of broadband UV radiometers at the “El Arenosillo” 
INTA station in Huelva (Spain) in September 2007 and June 
2011 [Vilaplana et al ., 2009]. These calibrations included 
spectral and angular characterization of the instruments and 
their absolute calibration, performed through the outdoor 
intercomparison with a reference Brewer spectroradiometer. 
Output voltages recorded by the UVB-1 radiometer were 
converted to UVER data applying conversion factors obtained 
from the “two-steps” calibration method [Seckmeyer et al ., 
1997; Webb et al ., 2006]. Anton et al. [201 la] compared the 
UVER data provided by the UVB-1 radiometer installed in 
Granada with radiative transfer model calculations for a 
cloudless sky; their results have shown high quality of the 
UVER data used in this paper. 

[s] The ground-based station is also equipped with a 
CM-1 1 pyranometer for measurements of global solar irradi- 
ance from 0.305 to 2.8 pm. This instrument is fully compliant 
with the highest ISO performance criteria with estimated rel- 
ative uncertainty better than 2% [Kratzenberg et al., 2006]. 
The stability of the pyranometer’s calibration has been peri- 
odically verified using a reference CM-1 1 instrument located 
at the study station and used only for inter-comparison 
purposes. The calibration factors showed variations below 


1% (four inter-comparisons performed between March 2005 
and June 2010) which guarantees the stability of the global 
solar irradiance data used in this work. 

[ 9 ] A Cimel CE-318 Sun photometer, co-located with the 
UVB-1 and CM-1 1 instruments, makes direct sun measure- 
ments with a 1.2° full field of view at seven wavelengths 
between 340 and 1020 nm, at every 15 min. This instrument 
is part of the Iberian network for Sun photometer aerosol 
measurements (RIMA), a scientific regional network feder- 
ated to NASA AERONET global network [Holben et al., 
1998]. From the direct sun measurements and using the cal- 
ibration constants provided by AERONET-RIMA, the aero- 
sol optical depths (AOD) at seven wavelengths are derived 
following the method described in the works of Alados- 
Arboledas et al. [2003, 2008]. Furthermore, the inversion 
procedure of Olmo et al. [2006, 2008] is utilized to retrieve 
columnar aerosol optical and microphysical properties such 
as single scattering albedo (SSA). This inversion code uses as 
input parameters AOD data derived from direct Sun pho- 
tometer measurements, and sky radiance measurements in 
the principal plane configuration. 

2.2. Satellite Data 

[ 10 ] The OMI satellite instrument is on board NASA EOS/ 
Aura platform launched in July 2004 [Schoeberl et al. , 2006]. 
This instrument consists of a nadir viewing push-broom 
spectrometer that measures solar backscattered radiation in 
the spectral range from 270 nm to 500 nm with a resolution of 
0.55 nm in the ultraviolet and 0.63 nm in the visible. The 
OMI instrument has a 2600 km wide viewing swath and it is 
capable of daily global contiguous mapping. The footprint 
size of satellite pixel is 13 by 24 km at nadir increasing up to 
~ 1 50 km off-nadir viewing directions. 

[ 11 ] The OMI surface UV algorithm (OMUVB) is based 
on the UV algorithm for Total Ozone Mapping Spectrometer 
(TOMS) instruments developed at NASA Goddard Space 
Flight Center (GSFC) [Krotkov et al., 1998, 2001]. This 
algorithm estimates surface UV irradiance from lookup tables 
(LUTs) obtained by a radiative transfer model using the 
OMI-derived total ozone, surface albedo and cloud infor- 
mation as input parameters [ Tanskanen et al., 2006, 2007]. 

[ 12 ] In this study, the OMUVB product used is OPEDRate 
(Overpass Erythemal Dose Rate). In addition, OMUVB data 
set contains the Lambertian Equivalent Reflectivity (LER) at 
360 nm which is used for cloud characterization. Addition- 
ally, we use the Aerosol Index (Al) calculated from 33 1 nm 
and 360 nm radiances which gives information about 
absorbing aerosols. All these OMI products are downloaded 
from the Aura Validation data Center site at http://avdc.gsfc. 
nasa.gov for the OMI station overpass data. 

3. Methodology 

[ 13 ] The inventory of Saharan desert dust events occurred 
at Granada from 2006 to 2010 is based on a published 
methodology using synthetic information from models, 
back-trajectories analysis, synoptic meteorological charts, 
satellite retrievals and surface data. Detailed information 
about this inventory can be found in the works of Valenzuela 
et al. [2012a, 2012b] who evaluated Saharan dust aerosol 
optical properties and its dependence on source region and 
transport pathways. 
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Figure 1. Atmospheric transmittance in the UV erythemal 
spectral range as a function of aerosol optical depth at 
380 nm in the slant path (AODs) for solar zenith angles 
smaller than 60°. The solid exponential curve represents the 
values given by equation (3). 

[ 14 ] In this work, the dust effect on the UVER is described 
in terms of the relative aerosol transmittance (TuverX rela- 
tive to no-aerosol clear sky UNER 0 [Krotkov et al . , 1998]: 

T UVER m 
UVER UVER 0 

In this expression UVER represents the erythemal mea- 
surements recorded during desert dust intrusions, and 
UVER 0 corresponds to the erythemal data for the same solar 
zenith angle (SZA) estimated from the empirical expression 
[ Mcidronich , 2007]: 

UVER 0 = a <*)‘g£) (2) 

where /*o is the cosine of the SZA and TOC is the total ozone 
column in Dobson Units (DU) provided by the OMI satellite 
instrument (OMT03 product using NASA TOMS V8 
retrieval algorithm [Bhartia and Wellemeyer , 2002]). This 
parameterized radiative model can be adjusted using avail- 
able experimental data, at the local site [Koepke et al . , 1998]. 
Anton et al. [201 lb] calculated the coefficients (a, b and c) in 
equation (2) for Granada site using UVER measurements 
during cleanest air conditions at the site. They validated this 
empirical model using measurements collected during a 
period not previously used for calculating fitting coeffi- 
cients. The results show a reliability of the empirical model 
(2), which estimates UVER 0 with a mean absolute bias less 
than 2.5%. 

[ 15 ] Finally, the atmospheric aerosol transmittance is cal- 
culated from equation (1) by using measured UVER values 
averaged within ±2 min of each Cimel AOD retrieval during 
desert dust events (2006-2010). 


[ 16 ] In order to analyze the effects of desert dust events on 
OMI-derived surface UV irradiances a single OMI ground 
pixel most closely collocated with Granada station is selected 
as the best match for each day. We used OMI pixels with 
centers from 1 km to 78 km from the study site, with the mean 
and median values being 17 and 11 km, respectively. Five 
surface UVER measurements within ±2 min from the OMI 
overpass at ~1 3:30 are averaged for comparison with the OMI 
data. Additionally, aerosol information derived from Cimel 
Sun photometer is averaged between 12:30 and 14:30 solar 
time on each day. 

4. Results and Discussion 

4.1. Effects on Ground-Based Measurements 

[ 17 ] The mean value (±1 standard deviation) of the atmo- 
spheric aerosol transmittance (equation (1)) is 0.89 ± 0.06 or 
1 1% average UVER reduction during desert dust intrusions 
over Granada. In ~12% dust cases UVER decreases more 
than 20% (Tuver < 0-8), which represents significant 
reduction in UVER due to aerosol absorption [Krotkov et al ., 
1998]. 

[is] Experimental studies about the evaluation of atmo- 
spheric transmittance in the UV range associated with desert 
dust episodes are very scarce: Diaz et al. [2007] reported 
atmospheric transmittance values between 0.92 and 0.95 
for global UV irradiance (280-363 nm) measurements in 
southeastern Spain and Kalashnikova et al. [2007] evaluated 
the atmospheric transmittance at 340-380 nm during 2 dust 
events in Australia, showing values larger than 0.90. These 
studies are focused on specific short-term dust events, while 
our study is based on a long-term inventory of dust episodes. 
Figure 1 shows measured UVER transmittance (Tuver) as a 
function of AOD at 380 nm in the slant path (AOD s ) for 
solar zenith angles (SZA) less than 60°. This variable takes 
into account the aerosol extinction in the slant column and 
it is derived by multiplying AOD with the air mass factor 
equal to sec(SZA) for each measurement [ Garcia et al, 
2006; Kazadzis et al., 2009]. From this figure, it can be 
seen that Tuver decreases as AOD s increases, with mini- 
mum Tuver values around 0.6 for AOD s of 1.5 (attenuation 
of ~40%) which highlights great influence of desert dust on 
measured UVER values. 

[19] The Tuver dependence on AOD can be parameter- 
ized as the exponential expression: 

Tuver = exp(-/r • AOD s ) (3) 

where k depends on aerosol single scattering co-albedo, \-u 
[Krotkov et al. , 1 998]. This parameter has been derived from 
the lineal regression analysis between the logarithm of 
Tuver and AODs, resulting in a mean value of 0.291 ± 
0.008 (R 2 ~ 0.7). Corresponding curve is shown in Figure 1 . 

[ 20 ] Krotkov et al. [1998] used a radiative transfer model 
to obtain the parameter k for different aerosol models. For 
non-absorbing aerosol (e.g., anthropogenic sulfate), they 
showed k values less than 0.15. Such small k values (i.e., 
Tuver ~ 1-2% for AOD = 0.1) are explained by the 
mutually compensating effects between the reduction of the 
direct flux by aerosol extinction and corresponding increase 
of the diffuse flux by aerosol scattering. For two dust 
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Figure 2. Correlation between OMI and ground-based 
UVER data for desert dust cases under cloud-free satellite 
pixels (in red). Subset of data for pristine (clear-sky) cases 
under cloud- and aerosol-free conditions are shown in gray. 
The study period is between January 2006 and December 
2010. The solid black line is the zero bias line, unit slope. 

models, the authors showed that the k values are substan- 
tially higher (between 0.3 and 0.6 for more absorbing dust). 
The larger k values can be explained by strong absorption of 
UV radiation by the mineral compounds of the desert dust 
such as hematite and goethite [Horvath, 1993; Alfaro et al. , 
2004]. The lower k ~ 0.3 value reported in our work com- 
pared to Krotkov et al. [1998] simulations can be explained 
by the location of our site thousand kilometers away from 
the desert dust source. Long-range transport of desert dust 
results in aging and mixing with other aerosol types leading 
to modification of its optical properties [Bauer et al . , 201 1]. 
Recently, Rodriguez et al. [2011] have shown that anthro- 
pogenic emissions from crude oil refineries and power 
plants, located in North African industrial areas, contribute 
to desert dust mixing with other type of particles. 

4.2. Dust Effects on Satellite Data 

[ 21 ] OMI surface UV algorithm assumes that clouds and 
aerosols are non-absorbing and, therefore, the satellite- 
derived surface UV irradiances are expected to show over- 
estimation for regions that are affected by absorbing aerosols 
such as smoke or desert dust [Tanskanen et al., 2007]. 

[ 22 ] The variability of cloudiness within the OMI pixel 
(13 by 24 km for nadir viewing) can lead to significant dif- 
ference between ground-based and satellite derived UVER 
data [Weihs et al., 2008; Anton et al., 2010]. To study the 
influence of desert dust aerosols on OMI UVER data, pre- 
dominantly cloud-free pixels should be considered in the 
analysis. We use OMI measured Lambertian Equivalent 
reflectivity (LER) values as proxy for cloud contaminated 
OMI scenes [Krotkov et al., 2001]. Specifically, we reject 
satellite scenes with LER values larger than 0.1 [Kalliskota 
et al., 2000]. Figure 2 shows the relationship between OMI 
retrieved and measured UVER data for desert dust events 


detected at Granada under cloud-free conditions (LER < 
0.1). The number of cloud-free days analyzed is 75 (69% of 
all dusty days). It can be seen that the correlation between 
satellite and ground-based UVER data is good (R 2 ~ 0.95), 
but a strong bias is also evident. The mean bias error (MBE) 
calculated as the average of the relative differences between 
OMI and measured UVER data (UVER omi -UVER exp / 
UVER exp ) is +22 ± 7% (OMI data being higher) where the 
uncertainty is characterized by the BE standard deviation. 

[23] In order to evaluate what part of this large OMI bias 
can be attributed to the presence of desert dust particles, the 
relationship between OMI and ground-based UVER data is 
analyzed for dust-free and cloud-free (i.e., pristine) condi- 
tions, which we call “clear sky.” Three different criteria are 
simultaneously applied for selecting clear sky cases. First, 
LER values smaller than 0.1 allow identify cloud-free sat- 
ellite pixels. Second, the clearness index (k t ) was used to 
characterize the atmospheric turbidity during satellite over- 
pass. The index is obtained from the ratio of the global solar 
irradiance to the extraterrestrial global solar irradiance on a 
horizontal surface [Alados-Arboledas et al., 2000]. In order 
to select cases with low turbidity, a conservative threshold of 
k t equal to 0.75 was chosen instead of the value of 0.65 used 
by other authors [e.g., Kudish et al., 1983; Udo, 2000]. This 
higher threshold guarantees that the selected cases corre- 
spond to the cleanest air conditions that occur at Granada. 
The third criteria was to select those cases with AOD at 
440 nm smaller than 0. 1 . Thus, we have implicitly assumed 
in this work that the atmospheric aerosol detected on clear- 
sky cases is the natural background. The pairs of satellite and 
ground-based UVER data recorded during clear-sky condi- 
tions are added to Figure 2. From this plot, it is highlighted 
that OMI bias is substantially reduced during satellite over- 
pass under clear-sky conditions. Thus, the MBE decreases to 
(+14.2 ± 4.1)% for these clear cases. On average, desert dust 
intrusions over Granada cause an increase of the OMI bias 
by additional 8 percentage points which is mainly related 
to the fact that current OMUVB algorithm assumes no 
absorbing aerosols [Tanskanen et al., 2007]. The main result 
of this assumption is the UVER overestimation due to the 
neglected aerosol absorption. In addition, since desert dust 
particles also reduce backscatter radiation reaching the sat- 
ellite, the OMI algorithm underestimates the effective cloud 
optical depth which produces an additional overestimation 
in UV radiation products [Krotkov et al., 1998, 2001]. This 
effective cloud optical depth is derived from matching the 
measured 360 nm radiance at the OMI overpass time with 
the modeled radiance assuming non-absorbing Cl cloud 
layer [Krotkov et al., 2001]. Nevertheless, we would like to 
point out that the analysis for clear-sky conditions shows a 
residual positive OMI bias around 14% which is not due to 
aerosol absorption, but can be related to several sources of 
uncertainty both in satellite and ground-based data. 

[ 24 ] Several papers [e.g., Krotkov et al. , 2005; Arola et al. , 
2005, 2009; lalongo et al., 2008; Kazadzis et al., 2009; 
Cachorro et al., 2010] have used the column aerosol 
absorption optical depth (AAOD) to quantify the error in the 
OMI UVER product due to the omission of correction for 
absorbing aerosols in the current OMI UV algorithm. This 
variable is calculated from the following expression: 

AAOD (A) = AOD( A) • [1 - S$4(A)], (4) 
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Figure 3. Ratio of OMI to ground-based UVER data as a 
function of the aerosol absorption optical depth (AAOD) 
extrapolated to 310 nm from Cimel values at 440 nm and 
870 nm using Absorbing Angstrom Exponent for desert dust 
cases detected between January 2006 and December 2008. 
The solid line represents the regression line, with the slope 
of the fit b = 2.69 (b parameter in equation (6)). 


where SSA is the single scattering albedo which can be 
retrieved from the sky radiance measured by the Cimel Sun 
photometer at 440, 675, 870 and 1020 nm [Olmo et al . , 
2006, 2008] or direct to diffuse irradiance ratio measured 
by UV-MFRSR [Krotkov et al ., 2005]. 

[25] For the OMI bias correction it is necessary to obtain 
the AAOD at UV wavelengths that are not currently mea- 
sured with Cimel data. In order to estimate AAOD at 3 1 0 nm, 
we use the following power law wavelength dependence 
[Bergstrom et al ., 2007]: 

AAOD( \) = c\~ AAE , (5) 

where AAE is the Absorption Angstrom Exponent which has 
been derived from the AAOD at 440 and 870 nm. 

[26] Figure 3 shows the ratios of OMI to ground-based 
UVER data plotted against AAOD at 310 nm estimated 
using equation (5). All dusty cloud-free data recorded for the 
period January 2006-December 2008 are included in this 
plot. It can be noticed that the ratio increases with increasing 
AAOD, confirming that dust aerosol absorption can partially 
explain positive OMI UVER bias found in the satellite- 
ground-based comparison. The regression analysis provides 
a slope of the fit of 2.1 ± 0.4, indicating the way in which 
the OMI bias increases with an increasing amount of aerosol 
absorption. Additionally, the linear least squares fit shows an 
intercept value of 1.14 ± 0.02 which corresponds to the 
remaining bias under pristine cloud- and aerosol-free con- 
ditions. This value suggests that OMI UVER data are biased 
14% high compared to the ground-based measurements for 
clear sky cases. 


[ 27 ] Based on the above results, OMI UVER data can be 
post-corrected using the expression proposed by Krotkov et al. 
[2005]: 


UVER OMI 

j TI/udOMI operational 

corr ~ (1 +b AAOD)' 


( 6 ) 


where the denominator accounts for the presence of absorbing 
aerosols during OMI overpass time, with b being the slope of 
the regression analysis performed in the previous paragraph. 

[ 28 ] Figure 4 shows the relationship between the reference 
ground-based UVER measurements and corrected OMI data 
for an independent data set corresponding to the period 
January 2009-December 2010 (not previously used for cal- 
culating the b parameter). It can be seen that correction 
method produces a clear reduction of the OMI bias. Thus, 
the MBE decreases from (+21 ± 5)% for operational satellite 
data to (+13 ± 4)% for corrected data, slightly smaller than 
the bias obtained for clear sky conditions (14%) which 
shows the level of improvement that may be reached with 
the off-line correction methodology. 

[ 29 ] For the off-line aerosol correction of the OMI UV 
data, Arola et al. [2009] used the equation (6) with b 
parameter equal to 3 and monthly AAOD values from the 
global aerosol climatology of Kinne [2009]. The improve- 
ment in the OMI UV that can be achieved with this correction 
has been also evaluated by comparing with experimental 
UVER data recorded at Granada for the period January 
2009-December 2010. The MBE (+13 ± 5)% is close to the 
bias obtained with the correction method used in our work, 
but with a larger standard deviation. Therefore, the Arola’s 
correction method could be successfully applied to post 
correct operational OMI UV products over geographical 



Figure 4. Correlation between OMI and ground-based 
UVER data for evaluation data set: desert dust cases 
detected between January 2009 and December 2010. Opera- 
tional satellite data are in red and corrected satellite data in 
gray. The solid black line is the zero bias line, unit slope. 
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regions where experimental aerosol measurements are not 
available. This correction is expected to be implemented in 
the next re-processing of the OMI data set. 

[ 30 ] The predecessor of the satellite OMI instrument was 
the TOMS instrument whose archived UV products were in 
fact corrected for aerosol absorption effects (dust and car- 
bonaceous particles) [Krotkov et al , 1998]. This correction 
was applied for those aerosol cases, identified from the 
TOMS measurements: a combination of the Aerosol Index 
(Al) larger than 0.5 and the LER less than 0.15. Thus, the Al 
= 0.5 is the threshold value chosen by TOMS UV algorithm 
to distinguish between absorbing and non-absorbing aerosol 
in the free troposphere. Contrary to the TOMS UV products, 
the current surface UV OMI algorithm does not use any 
correction for absorbing aerosols [ Tanskanen et al ., 2007]. 
For the 75 dusty cloud-free days selected in our work, the Al 
given by OMI is larger than 0.5 in 56 days (75%) which 
indicate that this satellite product may provide better infor- 
mation about the absorbing properties of the aerosol load in 
each pixel. However, the ratio satellite/ground-based UVER 
data presents a poor correlation with Al values (plot not 
shown) suggesting that Al itself is not the best quantity to 
evaluate the effect of absorbing aerosols on the OMI UV 
bias. It is better using AAOD (equation (4)) estimated from 
models [Arola et al. y 2009] or ground measurements 
[Krotkov et al . , 2005]. 

5. Conclusions 

[ 31 ] A long inventory of Sahara desert dust events recor- 
ded at Granada (Spain) has been used to analyze the influ- 
ence of this type of particles on broadband surface UV 
irradiance weighted by erythemal action spectrum (UVER) 
as retrieved by satellite (OMI) and measured by ground- 
based instruments. 

[ 32 ] The presence of desert dust aerosols over the study 
site causes average reductions of the UVER by about 1 1% 
with respect to clear-sky conditions. Reductions larger than 
20% are found in 12.5% of all desert cases. These results 
reveal that the desert dust particles markedly affect the 
propagation of the UV radiation through the atmosphere. 

[ 33 ] The UVER data derived from the OMI satellite 
instrument are biased high compared the ground-based 
UVER measurements during the desert dust cases with a 
mean relative difference of 22%. The analysis of pristine, 
clear-sky cases shows that 8% of the bias can be attributed to 
the fact that current OMI UV algorithm assumes no 
absorbing aerosols. Therefore, the effect of desert dust 
events on the UV irradiance derived from the OMI instru- 
ment cannot be neglected for regions like southern Spain, 
where the intrusions of the desert dust are frequent. 

[34] The aerosol absorption bias can be corrected off-line. 
The post-correction has been tested using an independent 
data set and resulted in reduction of the bias from ~21% for 
operational satellite UVER data to ~13% for corrected data. 
The remaining positive bias (OMI being higher), indicate 
additional sources for discrepancy. 

[ 35 ] Thus, OMI -derived UVER data have been shown to 
be overestimated in locations affected by desert dust. 
Therefore, reliable estimates of UV in these locations are 
dependent on the availability of quality assured ground- 
based measurements. 
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